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Abstract 

In this first part of a larger review undertaking the results of the first author and a 
part of the second author doctor dissertation are presented. Next we plan to give 
a survey of a nowadays situation in the area of investigation. Here we report on 
what follows. 

Calculation of the partition function for any vector Potts model is at first re- 
duced to the calculation of traces of products of the generalized Clifford algebra 
generators. The formula for such traces is derived. 

The latter enables one , in principle, to use an explicit calculation algorithm for 
partition functions also in other models for which the transfer matrix is an element 
from generalized Clifford algebra. 

The method - simple for Z% case - becomes complicated for Z n , n > 2, however 
everything is controlled due to knowledge of the corresponding algebra properties 
and those of generalized cosh function. 

Hence the work to gain the thermodynamics of the system, though possibly tanta- 
lous, looks now a reasonable, tangible task with help of computer symbolic calcu- 
lations. 

The discussion of the content of the in statu nascendi second part is to be found at 
the end of this Part 1 presentation. This constitutes the section V. 

PACS numbers: 05.50.+q, 

Keywords: Potts models, generalized Clifford Algebras 
I. Introduction 

The main idea of all calculations to follow (see [8]) is to consider the task of 
determining the complete partition function for noncritical Potts models - 
as a problem from the theory of generalized Clifford algebras classified 
in [15]. 
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As a matter of fact, the transfer matrix approach has led the authors of 
[1,14 except for 20,21] to these very algebras, although this observation does 
not seem to be realized by the mentioned authors. 

In general - the transfer matrix technique for a statistical system (or - a 
lattice field theory) with the most general translational invariant and glob- 
ally symmetric Hamiltonian (Action) on a two-dimensional lattice - does 
generate appropriate algebras of operators which are of the type of algebra 
extensions of some groups [8,18,24]. 

If one considers Z n cyclic groups as symmetry groups of Hamiltonian (chi- 
ral Action) then the algebra generated by transfer matrix approach is the 

(n) 

corresponding Cj, algebra, and it is due to its properties, that the given 
model has the duality property [18]. 

For duality property we refer the reader to the review [19], and for Potts 
models, in general - to [2,23]. 

Calculations to be carried out here for the Z n - vector Potts model simplify 
tremendously in the case of n = 2 i.e. for the Ising model and there lead 
to the known complete partition function [5] (see also [22] for modern pre- 
sentation) which after carrying out the thermodynamic limit, goes into the 
Onsager formula [16,17]. 

The method we choose to proceed with, is an appropriate generalization of 
the one used in [22] which consists there in reducing the problem of finding of 
the partition function for the Ising model to calculation of Tr (Pi...P s ), where 
Tr is the normalized trace while P's are linear combinations of 7 matrices - 
generators of usual Clifford algebra naturally assigned to the lattice. 
Then the observation that Tr {P\...P2 S ) is just a Pfaffian [3] of an antisym- 
metric matrix formed with scalar products of P's leads one to calculation of 
the determinant from this very matrix. 

The method proposed in [5,22] is purely algebraic, and though probably 
not the shortest one, it lacks ambiguities of other methods and uses well 
established, simple language of Clifford algebras. For other aspects and 
connotations of such an approach see [9,10,11,12]. 

Our paper is organized as follows: 

In the second section we write down the Z n -vector Potts model in a form 
resembling (and generalizing!) the Ising model without external field and 
then we represent the transfer matrix (an element of C^p !) as a sum of 
expressions proportional to Tr (7^. ..7^), where this time 7's are generalized 
7 matrices. 

Equivalently, the Hamiltonian for this Potts model can be looked upon as 
an Action for the Z n chiral model on the square lattice. 
In the third section we derive the formula for Tr (7i 1 ..-7j s ). 

(n) 

The last section is to supply an inevitable information on Cf, algebras 
for the Reader's convenience as well as some calculations avoided in the 
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current text in order to make the presentation more transparent. 



We would like to end this introduction by quotation from R. J. Baxter's 
book (see [2] p. 454): " The only hope that occurs to me is just as Onsager 
(1944) and Kaufman (1949) originally solved the zero-field Ising model by 
using the algebra of spinor operators, so there may be similar algebraic meth- 
ods for solving the eight-vertex and Potts models" . 

Our suggestion then is that these very algebras are just generalized Clifford 
algebras, and the presented paper is aimed to deliver arguments in favor of 
that point of view. 



II. The transfer matrix as a polynomial in 7's 

In the following the transfer matrix M for the Z n -vector Potts model is repre- 
sented in a form of "multi-sum" of expressions proportional to Tr (7^. ..7^). 

Let us assign to the set of states for this Z n vector Potts model on a 
p x q torus lattice (p rows, q columns), a set 



Def: S = {{s i:k ) = (p X q) ; s ijfc € Z n ] , 

where we have chosen a multiplicative realization for the cyclic group Z n = 
(^O/Lo an< ^ Sik denotes a matrix element of the (p x q) matrix. Here, 
naturally, uj denotes the primitive root of unity. 



The total energy E is then given by: 

p,q , p,q 



i,k=l i,k=l 

(2.1) 



while the partition function is defined to be: 

z= E e x p{-«} (2.2) 

(«i,fc)es 

One sees that the total energy of the system as represented by (2.1) is at 
the same time - a generalization of its own Z2 Ising case and - the Action of 
the corresponding chiral model (for connection between lattice gauge theory 
and spin systems see [6] and references therein). 

The partition function could be written in terms of transfer matrix and for 
that purpose we introduce the following notation: 
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Notation: 



Sk 



i=l 



(Si,k) = (si,S 2 ,...,S q ) 



( s l,k \ 
S2,k 



( S U \ 



5 2,fc 



(2.3) 



V 4* / 



With the notation (2.3) adopted, the partition function Z may be now 
rewritten in a form 

f q 9 ( \ 

(2.4) 



k=i 



s\ ,...,Sq K k — 1 

after the natural periodicity conditions have been imposed. 



Periodicity conditions: 

s q+1 = si, (4)i = (4) p+ i ; k = l,...,q, 



(2.5) 



where [x) i denotes the i-th component of x. 

The matrix ^ 1 is a p x p generalized Pauli matrix with matrix elements 
Si+ij , where i, j G Z p = {0, 1, ...,p — 1} and "+" is understood as the Z p 
group action. 

We introduce also the u\ generalized Pauli matrix, which is one of the 

(n) 

three a±, cr^ - playing the same role in representing C\ generalized Clif- 

(2) 

ford algebras as the usual ones in representing the ordinary C 2v Clifford 
algebras via well known tensor products of a matrices [4]. 
It is now obvious that Z may be represented as 



Z = Tr M q 



(2.6) 



as we have 



Z= M(s 1 ,s 2 )M(s 2 ,s 3 )...M(s q ,s 1 ) , 



si,...,s q 



where matrix elements of the transfer matrix M are given by: 

M{s,f) =exp{26Re(V-^ i s)}exp{2aRe(s*-s')} . (2.7) 

It is convenient to consider the matrix M as a product M = BA, where the 
corresponding matrix elements are identified as 



A (a", a*) =exp{2aRe(s"* •a')} 
and 

B (s, s") = exp {2b Re {s* ■ J2i «)} 8 (a, s" 



(2.8) 
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As all these A,B,M matrices are multiindexed it is obvious that they might 
be represented either as tensor products of (n x n) matrices (p times) or as 
(n p x n p ) matrices. 
It is not difficult then to see that 

A = ® p a i.e. (2.9) 

A is the p-th tensor power of the (n x n) matrix a, which has the form of a 
circulant matrix W [ai\: 

n—l 

a = (o/,j) = (exp {2a Re (a/ -7 )}) = £ Vi = W [<n] , (2.10) 

1=0 

where I,JeZ n = {0, 1, 2, n — 1} and 

\l = exp|2aRe (u 1 ^ } . (2.11) 

In order to see that A and B matrices are just some elements of 
we shall express them in terms of operators and ; k = 1, 2, i.e. 
matrices typical for tensor product representation of generalized Clifford 
algebras via generalized Pauli matrices (see (A. 3)). 

Def: 

X k = 7(g) ... <g) 1(g) cri ® J® ... ®I (p-terms), 
Z fc = 7(g) ... (g)7(g)c73 (8)7(8) ... (8)7 (p-terms), 

where <ti and 0-3 are situated on the k-th site, counting from the left hand 
side. 

The matrix A may be therefore now rewritten as a product of (n p x n p ) 
matrices 

p n—l 

A=l[w[X k } , where W [X k ] = ^ \X l k . (2.12) 
k=l 1=0 

Similarly, for the matrix B we derive: 

( P ~1 



B = exp \ {Z k X Z k + x + Z'^Zk) \ , (2.13) 
where Z p+ \ = Z\. 



k=i 



The formula (2.13) follows from the simple observation that matrix elements 
of 
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Z k ~ 1 Z k+ i + Z^-yZk (multiindexed by ~~s and s") give exactly In of the cor- 
responding term of (2.8) expression for B. 

The 5 function arises due to the fact that u 3 = (Si^juj 1 ) (see Appendix) 
and the exponentiation of matrix elements is possible because B is simply 
proportional to unit matrix. 

Once A and B have been represented as in (2.12) and (2.13) it is easy 
to express them in terms of generalized 7 matrices. Introducing then the 
tensor product representation (A. 3) we get: 



Z k 1 Z k+1 = H Sfc+i for odd n , 

and 

Xk = e^-sr ^ 

Z k l Z k+1 = £7™ Sfe+i fo r even n , 

where k = 1, 2, ...,p — 1 and £ 2 = u . 

The corresponding expression on the boundaries - read: 



(2.14) 
(2.15) 



Z- 1 Z 1 = ?77^ _1 7i ^ odd n , (2.16) 

and 

Z p 1 Z 1 = C 1 U7$~ 1, n f o r even n , (2.17) 

where 

uj-U = ® p oi . (2.18) 
For the proof of (2.14)-(2.17) use (A.6) and (A.7). 

From now on we shall proceed with formulas for n being odd, loosing nothing 
from generality of considerations while corresponding formulas for the case 
of n even are easily derivable from those for the odd case. This in mind we 

get 

A = f[ W KSr 1 ^] , (2.19) 
k=i 

B = ex P [ b E (^r Wi + i k +l%) I x exp {buj;- 1 ^ + bu-^-%} ■ 

(2.20) 

Our first goal is then achieved if one notes that 

U = f[ir 1 7k (2.21) 

fc=l 

i.e. the transfer matrix M is now expressed in terms of generalized 7 
matrices. 
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Before proceed, it is rather trivial and important to note that U n = 1, 
Z% = 1, X% = 1, with obvious implication of the same property for the n-th 
order polynomials in (2.19) and (2.20). 

Our second and the main goal of this section is to represent the transfer 
matrix in a form - reducing the Tr M q problem to calculation of Tr {ji 1 Ji 2 ---Ji s 
for some collections of 7's. 

(Note that for n = 2 the way to get the complete partition function is 
shorter as there, it is enough to reduce the Tr M q problem to calculation of 
Tr (P\P2...P S ) where P's are linear combinations of 7's. Hence the number of 
necessary summations is much, much smaller than in the case n > 2, where 
it is rather useless to try to represent A and B matrices in that convenient 
form). 

For that to do we shall deal now with the matrix B, to reveal its adequate 
for the purpose - structure. 

Let us start with observations following from U n = 1 property of U. Define 
matrices (see (A. 11)) to be 



i=0 

71—1 

i=0 



(2.22) 



and also 



B+ =exp{&^7^ 1 7 i} (2 23) 

B k =exp{6w fc 7™ l7 y p } , where k = 0, 1, n - 1 . 



n-1 



Then we have 

exp{6f/7r 1 7i} = ^5+y+ (2.24) 

and 

n-1 

expj^-sr^} = e w • ( 2 - 25 ) 

k=0 

Equations (2.24) and (2.25) follow from examining exp{Ux} via series ex- 
pansions modulo n, as in (A.8)-(A.10). 

The matrices have an important property (see (A. 11)) 

[Vt] n = V± (2.26) 

hence from (2.20), (2.24), (2.25) and the commutativity of matrix arguments 
of B 

(U, 7^ _1 7fe + i... etc.) it follows that: 

n-1 

B = E ffiVkBrvD (2-27) 

Lk=0 
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where 



B+ = exp \b tT V+i J> £ fe + , (2-28) 

and 



a=l 



B k = exp 1 6 £ 7^+1^ | ^ • (2-29) 

Expression (2.27) for 1? becomes still simpler due to the remarkable property 
of y fc 's: 

V k Vi = ; fc / / (2.30) 

where, for the moment, = V^, (see Appendix for proof). 
Hence 

n-l 
fc=0 

- all terms of the k-th summand - commuting. 

As for the commuting of various matrices involved in representing of the 
transfer matrix, note that 

[A, B]^0, [U, A] = and [V+, A] = [V~,A] = . 

All this is sufficient to write: 

M^Y^lBtB-AyiV+VrY (2.32) 

1=0 

where 1 < r < n — 1 and V q = V r as V n = V. 

It is not difficult to see how "r" arises. Namely "r" is the residual of the 
quotient {fzij- 

This n — 1 multiplicity in formula (2.32) makes the problem of thermo- 
dynamic limit more interesting and involved. 
If one assumes however that 

assumption: 

q = n + l(n-l) , (2.33) 
with I an arbitrary integer, then 

n-l 

M q = J2i B ? B r A } 9v i +v r ■ (2-34) 

1=0 

Note that there is no multiplicity in (2.32) for the Z2 case (Ising model) and 
also for that case V fc + = =V k , k = 0,1. 

The more: V^V^ = V k and (2.34) reduces its form considerably. 
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Comparison: 

In order to compare (2.34) with the known expression for M q in the case 
of Ising model [22] one should note that, for Z2 , (2.15) and (2.17) differ 
only by i and —i from (2.14) and (2.16) correspondingly, hence we have 



M q = (B-A) q V + + {B + A) q V- , (2.35) 



where 



B - = B o B o = ex P | 26i ( £ 7a7a+l - 7p7l) 
B+ = B+B^ = exp !^2bi ^ £^7 a '7a'+i - %7i 



and 



V+ = V = X - (1 + U) , V- = V 1 = X - (1 - U) . (2.35a) 

The formula (2.35) coincides then with the one known for Ising model [22] 
apart from the obvious (see (2.1)) and insignificant scaling of constants a 
and b by factor 2. 

Using the formula (2.35), the notion of Pfaffian and its relation to determi- 
nant - the author of [22] reobtained the complete partition function leading 
to the famous Onsager formula. 



Having the same goal in mind we are going at first to examine the expression 
(2.34) in order to see how (the polynomial in 7's !) M q is represented as a 
multisum of summands proportional to Tr (7^ ...7^). For that purpose we 
write: 

B+ = exp (bpu+) , (2.36) 

and 

B- k = exp (bp- 1 ^) , (2.37) 

where 
i.e. 

u+ = W V^-V , (2.38) 



u; = J°(nr% , (2-39) 

while 

n 2 -l 

p = p(n)=u)^- (jfe = 0,l,...,n-l). (2.40) 

Therefore both B k and B^ become n — 1 order polynomials in u k and u k 
correspondingly (see (A. 10)). 

One also shows easily that v£ and v~ defined below 
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exp{67a Sa+i} =exp{6pw+} , 



(2.41) 



exp{6 7 "; 1 1 7 Q }=e X p{6p- 1 t;-} , (2.42) 

do satisfy: 

(O n = (O n = l a = l,...,p-l. (2.43) 

hence both expressions (2.41) and (2.42) become n — 1 order polynomials in 
matrices v+ and v~ correspondingly. 

n=2: 

For n = 2 the further job is extremely facilitated due to the fact that A 
becomes then of the form 



a = n p k Q k 



(2.44) 



fc=i 



where P's and Q's are some known linear combinations of ordinary 7 ma- 
trices and similar holds for B+, B- matrices from (2.35). 
The matrices V+, V 7 - from (2.35) have also simple form and thus Tr M q 
disentangles for n = 2 to be the sum only four summands of the Pfaffian 
type i.e. Tr (P^..^). 

n>2: 

Unfortunately this disentanglement is no more possible for n > 2 as the 
polynomial W (see (2.12)) no more is representable uniquely as the product 
of linear combinations of 7's and neither is B. 
Hence the multi-sum becomes more complicated. 

Nevertheless it is obvious that Tr M q problem reduces to calculation of 
Tr (7i 1 7i 2 ...7i s ) for some collections of 7's. 

n=2: 

In the case of Ising model, the four arising Pfaffians contribute to the par- 
tition function to give [22]: 



Z = 2 vq ' x \ if [ch2a'ch26' - sh2a' cos f (21 + 1) - sh26' cos f (2k + 1) 



+ 



k,l=l 
PA 



+ n ch2a'ch26' - sh2a' cos | (21 + 1) - sh26' cos ^f- 



+ 



+ J! ch2a'ch26' - sh2a' cos *f - sh26' cos | (2k + 1) 
k,l=i L 
p,q r 



+ 



-a n 

k,l=l 



ch2a'ch26' - sh2a' cos M _ s h2b' cos ^ 

9 p 



(2.45) 
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where a denotes the sign of T — T c and a' = 2a, b' = 2b. Both the square 
root and the <7-sign have appeared here because of the use of Pf 2 = det 
relation. 

n>2: 

Again, for n > 2, as we shall see in the following section, although the 
generalization of the Pfaffian is possible to the case of arbitrary n, its relation 
to any valuable generalization of determinant does not to be valid as the 
arising signum like function no more is an epimorphism of Sk onto Z n , except 
for n = 2 of course. (Sk - the symmetric group of fc-elemental permutations). 
However, one may write M q as the polynomial in 7's and then use the general 
formula from the following section. 

This representation of M q in terms of 7^. ..7^ products is given in the Ap- 
pendix. 

As a result we have the following structure of the complete partition function 
for the Z n vector Potts models: 

.. n— 1 

M<1 = -2 E E G (fl)n(3;ji,j2) , (2-46) 

j'i,j2=0n 6 r m 

where G(...) are known functions of parameters a and b, (see the Ap- 
pendix: (A. 15)) and 

TrQ = to i or Tr = , 

where i = i (u,j 1 ,j 2 ^ G Z' n = {0,1, ...,n - 1}. 

The dependence of i on its indices is easy to be derived using the most 
general, appropriate formula for Tr (7u7i 2 ---7j s ) supplied by the next section. 



III. Trace formula for any element of C\J 

In this section the explicit formula for trace of any element of algebra 
is delivered. 

The very formula is crucial for getting the complete partition function for 
Potts models and hence (see [2] p. 454) for solving several major problems 
of statistical physics being unsolved till now since many years. 

In) 

The problem of explicit trace formula for M q £ C^ p , is decisive in 
calculation of Z function for those models on the lattice in which the transfer 

(n) 

matrix is an element of C\ p . 

We proceed now to derivation of the very formula. 

Note! By definition, in this section Tr map is normalized i.e. Tr / = 1. The 
derivation has the form of a sequence of lemmas. 
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Lemma 1. 

Let k 7^ n modn , k £ N ; then Tr (7u--7i fc ) = . 0" 

Proof: The same as for usual Clifford algebras. Use the matrix U defined 
by (2.18). 

Lemma 2. 

Tr (7z 1 7i 2 -"7i fc ) 7^ iff there exists permutation 5 € S^n , such that 

V(l) = V(2) = ••• = *<t(ti)> V(n+1) = ••• = V(2n)j •■■ > V(fen-n+l) = ••• = 

Proof: The proof follows from observation that due to (A.l) if no n-tuple of 
the same 7's exists then Tr(...) = 0. Other steps of the proof are reduced 
to this first one. 

It is therefore trivial to note, but important to realize, that: 
Lemma 3. 

Tr (7ii---7i fc ) = or / E Z n - the multiplicative cyclic group of n-th roots 
of unity. 

In Lemma 3 A; is again an arbitrary integer while in all preceding lemmas, 
and in the following, i\,i2, ■■■,ik run from 1 to number of generators of the 
given algebra. This number was chosen to be even, however note [15] that 
the "odd case" problem is reduced to this very one due to the properties of 
generalized Clifford algebra representations. 

The major problem now is to determine this value "0 or /" for arbitrary 
set of indices i± , 12 , . . . , ik ■ 

In order to do that define a signum like function K (unfortunately it is an 
epimorphism only for n = 2) - as follows: 

Def: 

K:S p ^Z n ; e„(i)&„(2)-&a( P ) = K (a) Q^-Qp , 

where 0's satisfy (A.l) except for the condition 7" = 1, which is now re- 
placed by ef = 1. 

This definition being adapted, it is now not very difficult to prove: 
Lemma 4. 
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T*(7«-7vJ =K{£)K{a) , for 

a ) V(l) = ••• = V(n) j ••• > V(pn-n+l) = — = V(pn) and 

b) ia(n) < h(2n) < ■■■ < h(pn) , 

where a = S o a , while S is a permutation of the elements {n, 2n, ...,pn}. 
(The group of S's is naturally identified with an appropriate subgroup of 
Spn)- 

Proof: The proof relies on observation that these are only different n- 
tuples which are "rigidly" shifted ones trough the others, i.e. there is no 
permutation within any given n-tuple. <0> 

The generalization of Lemma 4 to the arbitrary case of some of the n-tuples 
being equal - is straightforward. (The necessary change of conditions a) and 
b) is obvious). 

This in mind and from other lemmas we finally get: 



Theorem: 

0~(zS pn p SGSfi 

(*ct(piti+1); *ct([pi+P2]«)) X X ^ i}a(pn— p ; n+l) j •• ■■> *5(pn)) j 

with the notation to follow. <C> 
Notation: 

l 

V = (Pi,P2, -,Pi), Pi > 1, Pi = P, o = s ° u i and Sp is a subgroup of 

i=i 

S'pn isomorphic (for example!) to the group of all block matrices obtained 
via permutations of "block columns" of the matrix 



I P2 n 



where Ik is the (k x k) unit matrix. 



S - here denotes the multi-indexed Kronecker delta i.e. it assigns zero unless 
all its arguments are equal and in this very case S(...) = 1. The sum £' 
is meant to take into account only those permutations that do satisfy the 
conditions: 

a) a (1) < a (2) < ... < a (pin), ... , a (pn — pin + 1) < ... < a (pn), 
and 

b) a (1) < a (pin + 1) < ... < a (pn — pin + 1). 
Comments: 

1) For the case of n = 2 the theorem gives us the Pfaffian of the product 
7i 15 ...,7i 2 j as m the case, (and only! for n = 2) K (S) = 1 and we are left, 
as a result with only £' sum, while Kronecker deltas become functions of 
the same number of indices ij. 

2) The theorem solves our problem of TrM 9 , as any element of generalized 
Clifford algebra is a polynomial in 7's satisfying (A.l), including M q G . 
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IV. Final comments for the Part 1 of the presentation 

We have carried out our twentieth century investigation for the Z n vector 
Potts model known also under the name of planar Potts model. 
The similar investigation of the other Potts models, i.e. standard Potts 
models with two-site interaction [23] and multisite interactions as well, is 
being now carried out. 

However, it is to be noted here that the model considered in [14] possesses 
transfer matrix M = BA, where matrix A is a particular case of the one 
defined by (2.19) while B, though also expressed by matrices defined in 
section II, has a different polynomial (in these operators) in the expotential. 

As for the multisite interactions, the algebras to be used are the universal 
generalized Clifford algebras, introduced in [7]. Needless to say that these are 
standard Potts models which are of more interest because of their relation 
to a number of outstanding problems in lattice statistics [2,23]. 
To this end let us express our suggestion that the models of lattice statistics 
could be adapted (thanks to specific interpretation) to the domain of urban 
economics which, using the notion of entropy and information introduces as 
a matter of fact a kind of thermodynamics [13]. 

V. An outline of the second planned part content 

The content of the preceding chapters-except for the algorithm for calcula- 
tion of the complete partition function as presented above - was already pub- 
lished in twenty first century [25], [26]. This especially concerns the Clifford 
algebra technique omnipresent here and planned to play similar leading role 
in the next part of our review. The content of [26] from 2001 indicates the 
idea and a way how to use generalized Clifford algebra for chiral Potts models 
on the plane. The incessantly growing area of applications of Clifford alge- 
bras and naturalness of their use in formulating problems for direct calcula- 
tion entitles one to call them Clifford numbers. The generalized "universal" 
Clifford numbers are here introduced via k-ubic form replacing quadratic 
one in familiar construction of an appropriate ideal of tensor algebra. One 
of the epimorphic images of universal algebras k — C n = T(V)/I(Qk) is the 
algebra Cl n (k) with n generators and these are the algebras to be used here. 
Because generalized Clifford algebras Cl n (k) possess inherent Z^ grading - 
this property makes them an efficient apparatus to deal with spin lattice 
systems. This efficiency is illustrated in [26] by derivation of two major 
observations. Firstly, the partition functions for vector and planar Potts 
models and other model with Z n invariant Hamiltonian are polynomials in 
generalized hyperbolic functions of the n-th order. Secondly, the problem of 
algorithmic calculation of the partition function for any vector Potts model 
as treated here is reduced to the calculation of traces of products of the gen- 
erators of the generalized Clifford algebra. Finally the expression for such 
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traces for arbitrary collection of generator matrices is derived in [26]. 
Since the same 2001 year, due to the authors of [27] we know the form of the 
fc-state Potts model partition function (equivalent to the Tutte polynomial) 
for a lattice strip of fixed width and arbitrary length. 

From 2005 year Alan D. Sokal 54 pages review [28] aimed for mathemati- 
cians too one may learn that "' the multivariate Tutte polynomial (known 
to physicists as the Potts-model partition function) can be defined on an 
arbitrary finite graph G, or more generally on an arbitrary matroid M, and 
encodes much important combinatorial information about the graph'". 

Alan D. Sokal discusses there '"some questions concerning the complex 
zeros of the multivariate Tutte polynomial, along with their physical inter- 
pretations in statistical mechanics (in connection with the Yang-Lee ap- 
proach to phase transitions) and electrical circuit theory.'" Quite numerous 
open problems are also posed in [28]. For many references see both [27] and 
[28]. 

Coming back for a while to generalized Clifford algebra we mark their being 
just only mentioned in the abstract of Baxter paper [29] in which: '"The 
partition function of the N-state superintegrable chiral Potts model is ob- 
tained exactly and explicitly (if not completely rigorously) for a finite lattice 
with particular boundary conditions'". "' The associated Hamiltonian has 
a very simple form, suggesting that may be a more direct algebraic method 
(perhaps a generalized Clifford algebra) for obtaining its eigenvalues." ' Bax- 
ter - including his famous book [1982 Exactly Solved Models in Statistical 
Mechanics] comes back several time to the idea of Clifford or Clifford-like 
algebras' potential importance for the still not solved problem of complete 
partitions function obtained in a way Ising model was splved with help of 
Clifford algebras properties being used in a natural and elegant way. Man- 
ageable in a understandable way. Let us quote after Baxter from [30]:"' 
"This is rather intriguing - we are in much the same position with the chi- 
ral Potts model as we were in 1951, when Professor Yang entered the field 
of statistical mechanics by calculating MO for the Ising model. So on the 
occasion of his 70th birthday we are able to present him not only with this 
meeting in honor of his great contributions to theoretical physics, but also 
with an outstanding problem worthy of his mettle. Plus a change, c'est la 
meme chose." 'Here and there 

P. P. Martin's book [31] and his many subsequent papers in twentieth as 
well in twenty first century [32] are inevitable source of ideas and inspiration. 
And here now comes the 2005 year. 2005 was declared by some authors to 
be the major breakthrough for the chiral Potts model. This concerns 

an outstanding problem of the order parameters. See Baxter again [33] 
and [34] . There in [33] and then in [34] Baxter deals again with the problem 
of the order parameter in the chiral Potts model. He recalls that an elegant 
conjecture for this was made in 1983 and that it has since been successfully 
tested against series expansions, but as far as the author of [34] is aware 
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there is as yet no proof of the conjecture. 

2005. Again Professor Baxter. Here is an abstract of his Annual Confer- 
ence 2005 of the Australian Mathematical Society public talk .: entitled Lat- 
tice models in statistical mechanics: the chiral Potts model:" 'There 
are a few lattice models of interacting systems that can be solved exactly, in 
the sense that one can calculate the free energy in the thermodynamic limit 
of a large system. The interesting ones are mostly two-dimensional, such 
as the Ising model and the six and eight-vertex models. A comparitively 
recent addition to the list is the chiral Potts model. This is more difficult 
mathematically than its predecessors. While its free energy was calculated 
in 1988, until now there has only been a conjecture (a very elegant one) for 
the order parameters, i.e. the spontaneous magnetizations. This conjecture 
has now been verified, and in this talk I shall discuss the difficulties encoun- 
tered and the method used. "' The solution from 1988 he is referring to 
apparently refers to papers such as [35], [36], [37], [38] and others later, see 
[39]. By the way? - the authors of [39] implicitly indicate Generalized Clif- 
ford Algebras - in a footnote (3) referring there to Morris papers from 1967 
and 1968 [quoted in all Kwasniewski papers on subject]. The 69 pages, 30 
figures reviw [39] was written in honor of Onsager's ninetieth birthdayalso 
in order to present "'some exact results in the chiral Potts models and to 
translate these results into language more transparent to physicists'". 

This is more or less what the second part is planned to be about. By no 
means it sholud include review of numerous cotributions of Professor F.Y. 
Wu including not only Potts models [see: 

http://www.physics.neu.edu/Department/Vtwo/faculty/wu.../wupubupdated81803.htm 
but such fascinating papers as [40] of specifically personal interest of the au- 
thors |http://ii.uwb.edu.pl/akk/publl.htm|. Papers published in Advances 
in Applied Clifford Algebras such as [42] and the papers by the authors are 
to be included in the second part of this review too. 

Appendix 

1. generalized Clifford algebra is defined [15] to be generated by 

7i) ■•■)72p matrices satisfying: 

7i7j=^7j7i, i<3, l? = 1 i h3 = l,2,...,2p. (A.l) 

The very algebra has - up to equivalence - only one irreducible and faithful 
representation, and its generators can be represented as tensor products of 
generalized Pauli matrices: 

<J\ = (Si+ij) , 02 = (uj l 5 i+ ij) , o- 3 = (uj l 5ij) , (A2) 

where i,j G Z' n = {0, 1, ...,n — 1} - the additive cyclic group. 
One easily checks, that {<Tj}^ do satisfy (A.l) for n being odd. 



16 



Let I denotes since now the unit (n x n) matrix and let 



71 = 03 <g> i" ® I ® ... ® J ® I , 

72 = ci ® o"3 / <8) ... ® / ® / 



7 P = cri ® cri ® cri ® ... ® (7i ® cr 3 , . . 

71 = 0-2 ® /® I® ••• ® ^® I , 



72 = 0"! ® 0"2 ® ^ ® ••• ® I ® ^ 



7p = O"! ® (7l ® (7l ® . . . ® (7l ® (72 , 

then {7i,7j , i,j = 1, do satisfy (A.l) with w replaced by w -1 , hence 

(A. 3) are generators of the algebra isomorphic to (isomorphism is given 
by <?i <-> T3 in (A. 3)) (This very (A. 3) representation was chosen for technical 
reason - we get, for example, in calculations of section II, the matrix U 
without coefficients etc.). 
It is also to be noted that for n being odd 

<x 3 = <"V 2 . (A4) 

The case of n being even leads to similar representation with o\ un- 
changed but o"2 and 03 now equal to: 

02 = {csi+u) , ^ = ^rv 2 , (a.5) 

where £ is a primitive 2n-th root of unity such that 

(A. 3) then with these appropriate for case n = 2v generalized Pauli matrices, 

(n) 

provides us with the same type representation of C\ as the one for the case 
n = 2v + l. 

One then easily proves that 

<j™- l <j 1 =(Jz l for n = 2v + l K ' 

and 

(A.l) 



(T3 1 cr 2 = £ 1 WGi > 
cr^Vi = 1 for n = 2i/ . 



2. In this part of the Appendix we derive one useful formula, necessary for 
section II. 

Let x be any element of an associative, finite dimensional algebra with 
unity 1. 
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Then 

n-1 

exp {x} = Yl fi ( x ) , where 

/* 0*0 = £ ^E+iji , i = 0,..,n-l. 

We now express these /j's in terms expotentials. 
For that to do it is sufficient to note that 

fi(wx)=u i f i (x) , i = 0,l,...,n-l. {A.9) 

The (A.9) reveals the Z n symmetry properties of these generalized "cosh" 
functions and we get from this set of relations 

1 n—l 

f t ( x ) = - u~ H ex P {^ x } , * = 0,...,n-l. (A10) 

n fc=0 



3. For considerations of the section II we need the following 
Lemma 

Let U be as x above and in addition let U n = 1. Then V defined as 
follows 

n-1 

v = h E U% ' has the P r °P ert y : , A n \ 

v n = v . 

Proof: For the proof, just note that for some a j's 

n-1 



i=0 



and both sides of this identity equation must be symmetric in U l monomials. 
One concludes therefore that at = dj , i,j = 0, 1, n — l, hence - counting 
the number of all arising summands - one arrives at the conclusion of the 
Lemma. 



4. Here, the proof of the (2.30) formula follows. 
Let k = l + r,0<r<n — 1. Then 

1 n—l n—l 

W = ^ £ £ "~ rh K^) n+l2 - 

i 1= i 2 =0 

Introduce now new summation indices: i\ and i = i\ + ii- 



18 



Then we have 

n— 1 n— 1 



j 1= j=0 

because the summation over ii gives zero. 

The proof for other (+, — ), (— , +), (— , — ) cases is the same. 

5. In this part of the Appendix we derive the multisum structure of the 
complete partition function Z, for the Z n - vector Potts model with any n > 2. 
Since now on we shall use the following abbreviations: 

= Z' n <g> Z' n (g) ... ® Z' n , (k summands) , 

where Z' n = {0, 1, n — 1}, and 

^ , ^ 
to 2 = p [n) = p . 

Recall also , ; k G Z' n and i>+, u~ ; a = 1, ...,p — 1, defined in section 
II. Then we have, according to (A. 10): 

n-l , p— 1 n— 1 

# fe + = E /ii (bp) (u+) 1 n E/.(6P) w - 

h=0 a=l s=0 , a 1 9 \ 

n-l , p-1 n-l , 

5* = E /b (fe^ 1 ) (% ) n E /* M* . 

«2=0 q'=1 t=0 

In order to manage with abundance of indices we introduce a further, 
deliberate notation. 

Notation: 

l e r 2 , f,se r p _! i.e. 

£ = (^i , h) where ^,?2 £ Z' n etc. 
We also define: 

p— l p— l 

] 



F (b; k, L, S, f) = u Hh+h) p h-h fh {hp) fh ( 6p -i) n /si (6p) p -* x n f ptj ; 

V 7 i=l .7=1 

h (l,s,t) = [r P -^) h (jr\) h u (rVi)' n 1 feH-f 



p— 1 p— 1 

i 

and 



5 fc = B+B- . 
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Hence we may write 

B k = Y, F ( & 5 fc > L, S, f ) S (Z, S, f ) (A 13) 

L,S,f 

where 

Ler 2 , s,f€r p _i. 

It is important to recall now that: 
n-l 

M« = J] [S fcj A] 9 y+y- and [B fc , A] ^ . 

fc=0 

We introduce therefore again an appropriate notation. 
Notation: 

A (a; I) = fl A ir (a) w - *'" , and 
f (/) = fi (T^rf , while 

(see (2.21)) 

r=l 

(recall also (2.22)). 

This being adapted we may write: 

A= J>(a;/)f (/) (A14) 

Investigation of the multisum structure of M 9 matrix eventuates in rather 
transparent form of it, if one, (for the last time!) introduces an overall index, 
or rather "multiindex" for all reappearing summations. 

Notation: 

II = (^k, L±, Lg, Si, S g , Ti, T g , fi, Ig^j , 
G (ll) = ft F {b- k, L r , S r ,T r ) A (a; / r ) , 

All this together taken into account leads to 

Mq = ^2 E E G(n)n(n,ji,j2) (Ai5) 

jij2=0n 6 r m 
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where II G T rn ; m = 3pq + 1. 
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